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Abstract. We recalculate the strong lensing probability as a function of the 
image separation in TeVeS (tensor-vector-scalar) cosmology, which is a relativistic 
version of MOND (Modified Newtonian Dynamics). The lens is modeled by 
the Hernquist profile. We assume an open cosmology with Qf, = 0.04 and 
= 0.5 and three different kinds of interpolating functions. Two different galaxy 
stellar mass functions (GSMF) are adopted: PHJ (Panter-Heavens-Jimenez, 2004) 
determined from SDSS data release one and Fontana (Fontana et al., 2006) 
from GOODS-MUSIC catalog. We compare our results with both the predicted 
probabilities for lenses by Singular Isothermal Sphere (SIS) galaxy halos in 
LCDM (lambda cold dark matter) with Schechter-fit velocity function, and the 
observational results of the well defined combined sample of Cosmic Lens All- 
Sky Survey (CLASS) and Jodrell Bank/Very Large Array Astrometric Survey 
(JVAS). It turns out that the interpolating function /j(x') = x/{l + x) combined 
with Fontana GSMF matches the results from CLASS/ JVAS quite well. 
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1. Introduction 

The standard LCDM cosmology is very successful in explaining the cosmic microwave 
background (CMB, see, e.g., [78]), baryonic acoustic oscillation (BAO, see, e.g., 
[57]), gravitational lensing (see, e.g., [35]) and large scale structure (LSS) formation. 
However, LCDM faces some fundamental difficulties. From the observational point 
of view, the challenges to LCDM arise from smaller scales. For example, the theory 
cannot explain TuUy-Fisher law and the Freeman law [551 US]- The most difficult 
ones are the satellites problem and cusps problem. The most key problems are, of 
course, the unknown nature of Dark Matter (DM) and Dark Energy (DE). Before 
CDM particles are detected in the lab, science should remain open to the prospect 
that DM (and for the similar reasons, DE) phenomena may have some deep underlying 
reason in new physics. 

There are several proposals for resolving DM and DE problems by modifying 
Newtonian gravity or general relativity (GR) rather than resorting to some kinds of 
exotic matter or energy. MOND [35] was originally proposed to explain the observed 
asymptotically flat rotation curves of galaxies without DM, however, it was noticed 
that MOND can also explain TuUy-Fisher law and Freeman law [46l |47] . It is believed 
that MOND is successful at galactic scales [84l|83 (but see [37] for satellites problem). 
The challenges to MOND arise from clusters of galaxies [71], in which, some kind of 
dark matter, possibly some massive neutrinos with the mass of 2ev, is also needed 
to explain the dynamics of galaxies [3]. MOND and its relativistic version, TeVeS 
[6], are only concerned with DM, remain DE as it is. By adding a f{R) term in 
Einstein-Hilbert Lagrangian, where R is the Ricci scalar, the so called f{R) gravity 
theory can account for DE [I] [12l [13l [57l [58l [76l [III [82] . Another interesting theory 
is Modified Gravity (MOG) [50], it is a fully relativistic theory of gravitation that is 
derived from a relativistic action principle involving scalar, tensor and vector fields. 
MOG has been used successfully to account for galaxy cluster masses [9], the rotation 
curves of galaxies (similar to MOND) pjOj, velocity dispersions of satellite galaxies 
[STj . globular clusters [52] and Bullet Cluster [llj, all without resorting to DM. Most 
recently, MOG is used to investigate some cosmological observations (CMB, galaxy 
mass power spectrum and supernova), and it is found that MOG provides good fits 
to data without DM and DE ^ . 

Any modifications to traditional gravity theory must be tested with observational 
experiments. Gravitational lensing provides a powerful probe to test gravity theory 
[75l [85] . It is well known that, in standard cosmology (LCDM), when galaxies 
are modeled by a SIS and galaxy clusters are modeled by a Navarro-Frenk- White 
(NFW) profile, the predicted strong lensing probabilities can match the results of 
CLASS/JVAS quite well [I3[T7l[TS[Ti[20l[3l[33[Mllli[ia[ll[MllS[63l^ 
[60l[63ll59l[68l[Sl[74l[83l[86l[87]. 

This paper is devoted to explore the strong lensing statistics in TeVeS theory. As 
an alternative to LCDM cosmology, TeVeS cosmology has received much attention in 
the recent literature, in particular in the aspect of gravitational lensing [S] [231 [51], for 
reviews see [T] [T^- Before TeVeS, strong gravitational lensing in the MOND regime 
could only be manipulated by extrapolating non- relativistic dynamics [69.. . 54. , in which 
the deflection angle is only half the value in TeVeS [91 . In TeVeS theory, it is now 
established that, for galaxy clusters, both weak and strong lensing need extra DM to 
explain observations [3, ^301 [73] , possibly neutrinos with the mass of ~ 2ev, like 
the dynamics of galaxies. The situation is better for galaxies, as will be shown in 
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this paper. In our previous paper [22] , as a first try to calculate the strong lensing 
probability as a function of the image-separation A9 in TeVeS cosmology, we assumed 
a flat cosmology with fl}, = 1 — fi^ = 0.04 and the simplest interpolating function 
IJ,{x) ~ min(l,a;). In this paper, we assume an open cosmology with f2f, = 0.04 and 
fl\ = 0.5 and three different kinds of interpolating functions. As for mass function, in 
addition to the PHJ GSMF [66] used in our previous paper, we also adopt a redshift- 
dependent Fontana GSMF [31]. Further more, the amplification bias is calculated 
based on the total magnification of the outer two brighter images rather than the 
magnification of the second bright image of the three images as did in our previous 
work [221. 



2. TeVeS cosmology and deflection angle 

Gravitational lensing can be used to test TeVeS in two aspects. First, in the non- 
relativistic and spherical limit, TeVeS reduces to MOND. The deflection angle of the 
light ray passing through the lensing object can be calculated in MONDian regime 
(this will be discussed later). Second, the distances between the source, the lens 
and the observer are cosmological and thus depend on the geometry and evolution 
properties of the background universe. As argued by Bekenstein [6l[89], the scalar field 
(p, which is used to produce a MONDian gravitational acceleration in non-relativistic 
limit, contributes negligibly to Hubble expansion. According to the cosmological 
principle, the physical metric takes the usual Friedmann-Robertson- Walker (FRW) 
form in TeVeS [5], 

dT^ = -C^dt^ + a{tf[dx^ + /i(x)(d02 ^ ^-^^2 ^^^2)]^ 

where c is the speed of light, a{t) is the scale factor and 

r X-i/2sin(Xi/\) (if>0) 
fK{x)=lx {K = Q). (2) 

[ (-if sinh[(-if {K < 0) 

As in general relativity (GR), we define the cosmological parameters: 

where pb is the mean baryonic matter density in the universe at present time to 
(redshift z = 0), Pcrit(O) = 3H^/{8TrG) = 2.78 x lO^/i^^oMpc-^ is the present 
critical mass density, and Hq — lOO/ikms^^Mpc^^ is the Hubble constant. We choose 
a(to) = 1- Since dx = cdz/H{z), the proper distance from the observer to an object at 
redshift z is D^{z) = c /g^[(l + z)H{z)]~^dz, where the Hubble parameter at redshift 
z is (known as Friedmann's equation) 



H{z) = -= HoVW + + ^k{1 + + r!A. (4) 



a 
a 

The comoving distance from an object at redshift zi to an object at redshift Z2 is 

the corresponding angular diameter distance therefore is 

D{zi, Z2) = — /k[x(zi, Z2)]. (6) 
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In TeVeS, the lensing equation has the same form as in general relativity (GR), 
and for a spherically symmetric density profile |89| 

13 = 6 - -jr-a, a{b) = / ^ — dl, 7 

Ds Jo c^r dr 

where 6 = b/D^ and a{9) are the source position angle, image position angle and 
deflection angle, respectively; b is the impact parameter; Dl, Ds and Dls are the 
angular diameter distances from the observer to the lens, to the source and from 
the lens to the source, respectively; g{r) = d<I>(r)/dr is the actual gravitational 
acceleration [here $(r) is the spherical gravitational potential of the lensing galaxy 
and I is the light path]. It is well known that the stellar component of an elliptical 
galaxy can be well modeled by a Hcrnquist profile 

pin = - — I— r^, (8) 

27rr[r + rh) 
with the mass interior to r as 

where Mq = /q°° inr^ p{r)dr is the total mass and r^ is the scale length. The 
corresponding Newtonian acceleration is gAr(r) = GM{r)/r^ = GAIo/{r + r^j)^. 
According to MOND [551 [7TJ [72] , the actual acceleration g{r) is related to Newtonian 
acceleration by 

g{r)n{g{r)/ao) = gN{r), (10) 
where iJ,{x) is the interpolating function and has the properties 

{X, for a; <C 1 

1, for 1 

and gq = 1.2 x lO^^cms^^ is the critical acceleration below which gravitational law 
transits from Newtonian regime to MONDian regime. The concrete form of a /Lt(x) 
function should be determined by observational data (e.g., the rotation curves of 
spiral galaxies) and expected by a reasonable scalar field theory (e.g., TeVeS). The 
"standard" function one usually takes is fj,{x) = x/Vl + x^, which fits well to the 
rotation curves of most galaxies. Unfortunately, if the MOND effect is produced by a 
scalar field (such as TeVeS), the "standard" ^(x) function turns out to be multivalued 
[50] , On the other hand, a "simple" function fi{x) = x/{l + x) suggested by Famaey 
& Binney [55] fits observational data better than the "standard" function and is 
consistent with a scalar field relativistic extension of MOND [SD] [73] . 

In order to explore a broad class of modified gravity models, Zhao and Tian [52] 
proposed a parametrized modification function 



1 _ 9 



Ka/o-o) 9N 



1+ ( — 

9N 



kn 

ao ' 



(12) 



in which, MOND gravity corresponds to A; = 1/2. Substituting equation ([TOj) into 
equation ([T2|) with k = 1/2, we have 
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which can be easily solved to obtain the usual form of the function for MOND 

-2/n 



It is easy to verify that the "simple" and "standard" fi function are approximated 
with high accuracy by equation ([M]) with n ~ 3/2 and n — 3, respectively [92]- The 
requirement for a physical and monotonic fi function limits the parameter n to the 
range of 1.5 < n < 2.0. In this paper, we consider three cases: n = 1.5, 2.0 and 3.0. 

Since the MONDian gravitational acceleration g is explicitly expressed in terms of 
the Newtonian acceleration g^, it is very convenient to use equation (|12p to calculate 
the deflection angle 

Jo r 

^^[1 + {^r/^Y'-Ai (15) 

& r (r + Th) gN 



By using r = b^y^ + (l/b)^ and 9 = b/Di^, we have 

aiO) . 0".207/.- f ^) ^ r ^ [1 + K/^N)"/T/" 

where M = Mq/M^ and = 7.64 x IO^'^/i^^Mq is the characteristic mass of galaxies 
[55] . and 

- - 2.38 f-^V ^ f yrT^ + 0.05^^V . (17) 
5N \c/HoJ M \ Di^ e J 

In equations and ITTl) . the image position angle 6* and the scale length rh are in 
units of arcsecond ( ) and Kpc, respectively. 

We need a relationship between the scale length rh and the mass M, which could 
be determined by observational data. First, the scale length is related to the effective 
(or half-light) radius Re of a luminous galaxy by = i?e/1.8 It has long been 
recognized that there exists a correlation between and the mean surface brightness 
(/e) interior to i?e [26j: -Re oc (/) jO-83±o.08 ^ ^^^^^ ^^le luminosity interior to Re (half- 
light) is Le = L/2 = 7r(/)ei?g, one immediately finds Re oc L^ '^^. Second, we need 
to know the mass-to-light ratio T = M/L (x LP for elliptical galaxies. The observed 
data gives p = 0.35 [5T]; according to MOND, however, we should find p « [72]. In 
any case we have 

LcxMi/(i+P). (18) 

Therefore, the scale length should be related to the stellar mass of a galaxy by rh — 
j[]\,fi-26/{i+p) ^ and the coefficient A should be further determined by observational 
data. Without a well defined sample at our disposal, we use the galaxy lenses which 
have an observed effective radius Re (and thus rh) in the CASTLES survey [SS], which 
are listed in table 2 of ^89]. The fitted formulae for rh are 

1.26 

Kpc, for p=0.0. 







1 0.72 












[l.24 









'^h = .... 1.26/1.35 ■ (19) 

Kpc, for p=0.35 

In later calculations, except indicated, we use the fitted formula of rh for p = as 
required by MOND. 
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3. Galaxy stellar mass function 

In LCDM cosmology, mass function of virialized CDM halos can be obtained in two 
independent ways. One is via the generalized Press-Schechter (PS) theory, the other 
is via Schechter luminosity function. In TeVeS, however, the PS-like theory does not 
exist. Fortunately, the stellar mass function of galaxies is available in the literature, 
including the one constrained by the most recent data [311 [66] . 

Before giving the galaxy stellar mass functions (GSMF) appeared in the most 
recent literature, it is helpful to derive a GSMF directly from the Schechter luminosity 
function and mass-to-light ratio. The Schechter luminosity function is 

,W.r(A)°„p(-A)^i, ,20, 

For i/L* = {M/A'Q^/^^+P^ implied by equation ([18]), wc have a GSMF 



s " + 1 1 



1 n 



m7 



dM 

m7' 



(21) 



While the average number density of galaxies , the slope at low-mass end a and the 
slope of mass-to-light ratio p may be easily found from the published observational 
data or assumptions, the characteristic stellar mass of galaxies M^, can be derived 
from 

/>oo 

Plum = ni,^Pcrit(0) = / M(l>{M)dM, (22) 



Jo 

where pium is the luminous baryonic matter density (note that pium ^ Pb)- The 
characteristic mass Mi, is 

^lumPcrit(O) , s 

^^*-0J(a+p + 2)- ^2^^ 

For example, for ((/)^, a, fJiu,„,p) = (0.014/i3Mpc"^, -1.1, 0.003, 0.35) from [41 , Af^ = 
6.56 X 10^"h~^MQ; for the same parameters except that p = 0.0 (MOND), Af* = 
5.56 X 10^°h-^MQ. 

Fortunately the parameters in equation (|2ip have been determined by recent 
observational data. By determining non-parametrically the stellar mass functions of 
96545 galaxies from the Sloan Digital Sky Survey Data (SDSS) release one, Panter, 
Heavens and Jimenez [BB] (PHJ, hereafter) give the GSMF [55] 

f MY f M\dM 
0(A/)dAf = 0.(-j exp (24) 

where, we use </'(Af) to denote the comoving number density of galaxies with mass 
between M and M -\- dAf, and 

(/)^ = (7.8 ±0.1) X 10"3/^3Mpc"^ 

a = -1.159 ± 0.008, (25) 



Mi, = (7.64 ± 0.09) X 10^°/i~2a/, 



0- 



Most recently, in order to study the assembly of massive galaxies in the high redshift 
Universe, Fontana et al. [5T| fFontana. hereafter) used the GOODS- MUSIC catalog 
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to measure the evolution of the GSMF and of the resulting stellar mass density up to 
redshift z ^ i. The GSMF they obtained is 

1 a(z) 

(j){M,z)dM = 



cxp 



M 



M^z) 



dM 



M^izY 



(26) 



where 



04z)=n*(l + z)< 



n. 







0.0035, n\ ^ -2.20 ±0.18, 



a{z) = do + diz, do = —1.18, di = 
M4z) = 10*=^°+*^i^+^^^*^'/i-2Mq, 



-0.082 ±0.033, 



(27) 



M* = 11.16, il/* = 0.17 ± 0.05, = -0.07 ±0.01 

It would be interesting to compare PHJ and Fontana GSMFs to the mass function 
of galaxies in LCDM cosmology when the galactic halos are modeled by SIS. The 
comoving number density of galactic halos with velocity dispersion between v and 
t; + di; m [22] is 



{v)dv = 



exp 



(28) 
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Figure 1. Comoving number density for PHJ (solid), Fontana (dotted) and SIS 
halos (dash). Since Fontana mass function depends on redshift, four cases with 
redshift z = 0.0, 0.5, 1.0, 2.0 are displayed. For comparison, wo normalize the three 
mass functions to the same value of characteristic mass 
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For comparison, we need to transform equation 
mass M 

rr2oo soOtt 



M = 477 



Psis{r)r dr 



from vefocity dispersion to hafo 

^200Pcrit(2:), (29) 







where r2oo is the virial radius of a galactic halo within which the average mass density 
is 200 times the critical density of the Universe Pcntiz). Substituting the well known 
expression psis('') = w^/(27rGr^) into equation ([29|) . it is easy to find 

M{z) = 6.58 X 10^ (k^i^)^ + + ^^^^ + flAr^^^h-^AlQ, (30) 

where fim is the matter density parameter (including dark and baryonic components) 
[42] , Equation (|30p means that at any redshift z we should have AI cx u'^, or for our 
purpose, another form 

M _ ( 

we thus have the galaxy mass function for SIS halos 



(31) 



m) - ^ (- 



(3-2)/3 



exp 



M 

m7 



/3/3- 



dM 

m7' 



(32) 



We plot PHJ and Fontana GSMFs in figure 1 together with the galaxy mass function 
for SIS halos (comoving number density). For SIS halos, we use (0*,a,/3) = 
(0.0064/^3Mpc"^-1.0,4.0) [M]. For comparison, we normalize the three mass 
functions to the same value of characteristic mass A/^ ~ 7.64 x IG^^/i^^A/q. Note 
that, for Fontana GSMF, the comoving number density of galaxies decreases with 
increasing redshift, as expected 



4. lensing probability 



Usually, lensing cross section defined in the lens plane with image separations larger 
than A6 is cr(> AO) = nDll3^^Q[Ae{M) ~ AO], where e(a;) is the Heaviside step 
function and /3cr is the caustic radius within which sources are multiply imaged. This 
is true only when A6{M) is approximately constant within /3cr, and the effect of the 
flux density ratio gr between the outer two brighter and fainter images can be ignored. 
Generally this is not true, readers are referred to [22] for details. We introduce a 
source position quantity (3q^ determined by 



9r 



0{f3) de{(i) 



(33) 



h<0<b 



13 dp 

where 6*0 = ^(0) < 0, the absolute value of which is the Einstein radius, and O^-r is 
determined by c?/3/d6' = for 6* < 0. Equation ([55)1 means that when /3g_, < /3 < /3cr, 
the flux density ratio would be larger than q-^, which is the upper limit of a well defined 
sample. Therefore, the source position should be within [3q^ according to the sample 
selection criterion. For example, in the CLASS/JVAS sample, < 10. 

The amplification bias should be considered in lensing probability calculations. 
For the source QSOs having a power-law flux distribution with slope 7 (= 2.1 in 
the CLASS/JVAS survey), the amplification bias is B{P) — jT'^^ [65,, where, in this 



paper. 



0_de_ 

Pdp 



h<e<e 



\ P dp j,^. 



(34) 
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is the total magnification of the outer two brighter images. In our previous work |22j . 
however, the amphfication bias is calculated based on the magnification of the second 
bright image of the three images. 

Therefore, the lensing cross section with image-separation larger than AO and flux 
density ratio less than and combined with the amplification bias B(J3) is [75l[T9l[22] 

(t(> A9, < q,) = 2ttDI X 

3,, 

for AO < Aeo, 

ri''-\l3)dl3, for A9o < A9 < A6l,,, ^^^^ 
for AO > A9q^ , 

where is the source position at which a lens produces the image separation A9, 
AOq — A9{0) is the separation of the two images which are just on the Einstein ring, 
and A9q^ = A9{Pq^) is the upper- limit of the separation above which the flux ratio of 
the two images will be greater than qr- 

The lensing probability with image separation larger than AO and flux density 
ratio less than q,-, in TeVeS cosmology, for the source QSOs at mean redshift Zs = 1.27 
lensed by foreground elliptical stellar galaxies is [19l [20l [211 HI] 

P{> A9, <qr)^ / / U z / 0(M, z)(l + z)V(> A9, < q,)dM, (36) 
Jo "z Jq 

We plot in Figure [2] and Figure [3] the numerical results of the lensing probability 
according to equation p6|l . In TeVeS (solid lines), we assume an open cosmology 
with ilb = 0.04 and = 0.5, as implied by fitting to a high-z Type la supernova 
luminosity modulus [89]. The lensing galaxy is modeled by Hernquist profile with 
length scale rh — 0.72(M/Mi,)^-^^Kpc for constant mass-to-light ratio as required by 
MOND [see equation ((19])]. The interpolating functions with three cases n — 1.5,2.0 
and 3.0 are considered (top-down) according to equation (|13p . In order to investigate 
the effects of MOND on strong lensing, we also calculated the probabilities (dotted 
lines) with no modification to gravitation theory (i.e., in GR) and without dark matter 
(i.e., lensing galaxy is modeled by Hernquist profile). In this case, two types of the 
fitted formulae for the length scale rh with p = and 0.35 (top-down) are adopted. In 
TeVeS and GR (with no dark matter), we adopt the GSMF as the mass function (mf), 
with mf=PIIJ in Figure [2] and mf=Fontana in Figure [3j As did in our previous work 
[22] , We recalculate the lensing probability with image separation larger than A9 and 
flux density ratio less than qr, in flat LCDM cosmology {flm = 0.3 and = 0.7), for 
the source QSOs at mean redshift Zg = 1.27 lensed by foreground SIS modeled galaxy 
halos [a [45] [49]: 

dDP(z] 1"°° 

Psisi> A9,<q,)^ dz / ' dvn{v,z)asis{v,z)B, (37) 

Jo Jv^g 

where n{v,z) = (j)(v){l + z)^, which is related to the comoving number density (f>{v) 
given by equation ([^5)1 . is the physical number density of galaxy halos at redshift z 
with velocity dispersion between v and v + dv [49j , 

asis(.,.) = 16.^nY^V (38) 




Figure 2. Predicted lens probability with an image separation angle > A9 and 
the flux ratio < Qr = 10. For TeVeS (solid line) and GR (no CDM and without 
modification of gravity, dotted line), we assume an open cosmology with Of, =0.04 
and Oa = 0.5, model the lens as the Hernquist profile and adopt PHJ GSMF (I24I I: 
for standard LCDM (dashed line), we assume a flat cosmology with = 0.3 and 
SIa = 0.7, model the lens as the SIS and adopt the mass function II28I I. For GR, 
we consider two different mass-to-light ratio types and thus the expressions of 
rji , see equation 11911 . For comparison, the survey results of CLASS / J VAS (thick 
histogram) are also shown. 



is the lensing cross section, 




(39) 

is the minimum velocity for lenses to produce image separation > A9 and B is the 
amplification bias. We adopt a, /3) = (0.0064/i3Mpc~^, 198kms"\ -1.0, 4.0) 

for early- type galaxies from [M]. A subset of 8958 sources from the combined 
JVAS/CLASS survey form a well-defined statistical sample containing 13 multiply 
imaged sources (lens systems) suitable for analysis of the lens statistics [Sni ISl EZl [3S] ■ 
The observed lensing probabilities can be easily calculated [TSllinilll] by Pobs(> ^6*) = 
A^(> A6')/8958, where iV(> A9) is the number of lenses with separation greater than 
in 13 lenses. For comparison, the observational probability Pobs(> A.6) for the 
survey results of CLASS/ JVAS is also shown (thick histogram). It would be helpful for 
us to figure out differences among models to summarize the values of the probabilities 
P(> A0 = 0.3") in the Table 1. 




Figure 3. same as Figure 2, except that GSMF=Fontana for TeVeS and GR. 



5. Discussion and conclusions 

We have calculated the lensing probability with image separation larger than a given 
value in an open, TeVeS cosmology. The results are sensitive to the interpolating 
function fi{x) and mass function (j){M,z). For a given GSMF (PHJ in Figure [2] and 
Fontana in Figure[3]), the lensing probability decreases with increasing value of n [given 
in equation (fTl]) ]. Obviously, for PHJ GSMF (Figure [2|), the lensing probabilities 
calculated in TeVeS (solid lines for three cases of interpolating functions) are too 
large at small lensing image separations compared with the results of CLASS/ JVAS. 
This unreasonable result is further confirmed, when we note that, even the lensing 
probabilities in GR cosmology (with no DM, dotted lines) are much larger than that 
in LCDM cosmology (dashed line) at small image separations. Actually, however, 
this result can be easly explained: at small mass-end (corresponding to small image 
separation), the comoving number density for PHJ mf is much larger than that 
for SIS halos (Figure [T]), which results in the corresponding lensing probabilities 
according to equation ([55]) . This is why in our previous work [55], we calculated the 
amplification bias based on the magnification of the second bright image rather than 
the total magnification of the two images considered. According to the resolution 
of CLASS/ JVAS, however, it is difficult to resolve the two images for small image 
separations. Therefore, in this paper, we calculate the amplification bias based on the 
total magnification of the outer two brighter images, as usually done in the literature. 
On the other hand, if we adopt another most recent mass function, Fontana 
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Table 1. The predicted values of the lensing probabilities P(> A9 = 0.3") for 
all models and Pobs(> AS = 0.3") for CLASS/JVAS. 



TeVeS 



GR 



LCDM CLASS/JVAS 



PHJ 

Fontana 



11=1.5 
-2.460 
-2.873 



n=2.0 
-2.513 
-2.926 



n=3.0 
-2.552 
-2.965 



p=0.0 
-2.601 
-3.021 



p=0.35 

-2.691 

-3.112 



-2. 



-2.838 



GSMF (Figure [3]), we find that the predicted lensing probabihties in an open TeVeS 
cosmology with the "simple" interpolating function fi{x) — x/{l + x) [i.e., n = 3/2 
in equation ([H])] match the observational data of CLASS/JVAS quite well. Similarly, 
this is reasonable when we note that the comoving number density of galaxies for 
PHJ GSMF is much higher than that for Fontana GSMF at small mass-end (Figure 
1). Clearly, the "standard" interpolating function fj,{x) — x/y/l + x^ [i.e., n = 3 in 
equation (|14p ] is ruled out due to its too low lensing rates at small image separations. 
Interestingly, this conclusion is in agreement with the most recent result of Sanders 
and Noordermeer [75] , who constrained the interpolating function with the rotation 
curves of early-type disc galaxies. 

In our calculations for deflection angle in TeVeS cosmology, we have fixed the value 
of the critical acceleration uq and modeled the lensing galaxies with the Hernquist 
profile, and the only free choice is the interpolating function fi{x). We note that the 
PHJ GSMF includes all types of galaxies, whereas the mf for SIS halos includes only 
early-type galaxies, this can partly explain the relativey low abundance of SIS halos 
compared with PHJ GSMF in figure 1. On the other hand, Fontana GSMF, like PHJ 
GSMF, also includes all types of galaxies, but its value is close to (when z = 0) or 
lower (for high z) than the mf for SIS halos. Therefore, the major uncertainty for 
lensing probability arises from the GSMF, which is independent of any gravitational 
theory and should be determined by observational data. Can we conclude from Figure 
2 and Figure 3 that Fontana GSMF is preferred and PHJ GSMF is ruled out? Recall 
that, in LCDM cosmology, there are also uncertainties for the mass function derived 
from the luminosity function, the equation (pS)) [TBI fiSl IMj . Actually, the parameters 
we adopted in equation (</>*, u^, a, /3) = (0.0064/i3Mpc"^, 198kms"\ -1.0, 4.0) 

[2], are selected so that the predicted lensing probabilities Psis(> 0.3") can exactly 
match the observed value Pobs(> 0.3"), i.e., Psis(> Afl = 0.3") = Pobs(> M = 

0. 3"), see Table 1. The most recent parameters derived from SDSS DR3 [24] is 
((/)^,'i;*,a,^) = (0.008/i3Mpc"^,161kms"\ 2.32, 2.67), however, this wiU not affect 
our results. One can see clearly from Figure 2 that, at larger image separations, 
the predicted lensing probabilities (dashed line) are well bellow the observed values, 

1. e., Psis(> Afl > 0.3") < Pobs(> AO > 0.3"). As a matter of fact, in strong 
lensing statistics, one usually compares the predicted cumulate lensing probability 
at the image separation of A9 = 0.3", and regards the under-estimates at larger 
image-separations to be unimportant. We note, however, that there is an infiexion at 
A9 = 1.16" for Pobs(> A^) calculated from the well-defined sample of CLASS/JVAS 
(thick histogram) , and there are no physical interpretations for the flat part of the line 
when AO < 1.16". Although the sample of CLASS/JVAS is "well-defined" , this should 
not include each detail such as the inflexion, and other observations, like [33] , would 
provide more information at A9 < 1.16". So we can reasonably guess that the correct 
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observational data should avoid the inflexion, and the trend of the rising probabilities 
with smaller and smaller image-separations should continue inward at A6' < 1.16". In 
this sense, the predicted lensing probabilities in TeVeS cosmology, with a PHJ GSMF 
and the "standard" interpolating function, match the observational data quite well as 
shown in Figure 2. 

We also note that, in Figure 3, the lensing probabilities in GR cosmology (with 
no DM, dotted lines) are much lower than the observational data. This imply that, 
as an alternative to CDM, MOND can sufhciently account for the strong lensing 
observations. 
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